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According to the Landau criterion for superfluidity, a Bose-Einstein condensate flowing with a 
group velocity smaller than the sound velocity is energetically stable to the presence of perturbing 
potentials. We found that this is strictly correct only for vanishingly small perturbations. The 
superfluid critical velocity strongly depends on the strength and shape of the defect. We quantita- 
tively study, both numerically and with an approximate analytical model, the dynamical response 
of a one-dimensional condensate flowing against an istantaneously raised spatially periodic defect. 
We found that the critical velocity v c decreases by incresing the strength of the defect Vo, up to to 
a critical value of the defect intensity where the critical velocity vanishes. 



Introduction. Dissipationless flow in homogeneous 
bosonic systems is one of the manifestations of the inter- 
play between interatomic interaction and Bose-Einstein 
condensation. Bogoliubov ground-breaking calculation of 
the many-body properties of a dilute condensate showed 
that this interplay produces a linear dispersion of the 
low-energy excitation modes which, according to Lan- 
dau, inhibits dissipation when the group velocity of the 
fluid is smaller than the sound velocity Q, 0, 01 • How- 
ever, the question if a Bose-Einstein condensed dilute gas 
(BEC) really displays this basic characteristics of super- 
fluidity remained open until its experimental realization 
in 1995. Eventually, several typical superfluidityphe- 
nomena like vortices with quantized circulation [j, |5( , 
the scissor mode 0, Q and Josephson dynamics p| have 
been experimentally observed. Instead, the direct exper- 
imental verification of the Landau critical velocity has 
been more elusive. The MIT group studied the response 
of a trapped elongated condensate perturbed by a blue 
detuned laser beam [j| . They found a transition between 
the superfluid and the dissipative regime at a laser ve- 
locity smaller than the estimated sound velocity. The 
reason for this discrepancy is not fully understood, but 
it might be related to the topological creation of vortices 
[Tfj|. It is worth noticing that even in superfluid 4 He 
the observed superfluid critical velocity is smaller than 
Landau's prediction [Tl|. 

In this paper we study the response of a one dimen- 
sional homogeneous condensate flowing with momentum 
p against an instantaneously raised defect. The analysis 
is done numerically, solving the Gross-Pitaevskii equa- 
tion (GPE), and analytically, with both a dynamical Bo- 
goliubov approach and a nonlinear two-mode approxima- 
tion. The Bogoliubov linear framework provides a clear 
dynamical analysis of the Landau criterion, predicting 
the transition from a superfluid to a dissipative regime 
at a condensate group velocity equal to the sound veloc- 
ity. However, and this is the main result of this paper, we 
show that nonlinear corrections to the Landau criterion 
become very important even for small strengths of the 
defect. 



We consider, as a defect, a spatially periodic potential. 
This choice is convenient for two reasons: first, in the 
linear regime it is possible to selectively excite a given 
pair of quasi-particle modes, and study their behaviour 
below and above the critical velocity; second, periodic po- 
tentials can be easily taylored and accurately controlled 
exp er iment ally. 

Our analysis has several analogies with the Bragg 
spectroscopy of condensates, which has been exten- 
sively invest igat ed both theoretically and experimentally 
E E Q CI 03- Thc crucial difference is that we 
focus on the signature of the energetic instability associ- 
ated with the Landau criterion of superfluidity, and on 
the study of how the nonlinear dynamics modifies the 
Landau critical velocity. 

In particular, for a spatially periodic defect, we find 
that the Landau critical velocity v c rapidly decreases with 
increasing defect strength, up to a critical value where 
superfluidity disappears. This result might be related to 
the downwards shift of the excitations frequencies with 
respect to the Bogoliubov spectrum, as predicted in |15| . 
While the regimes investigated in Bragg spectroscopy ex- 
periments are still described with a good accuracy by 
linear Bogoliubov theory, the experimental investigation 
of our predictions would require Bragg experiments in 
highly non-linear regime, i.e. beyond the regimes inves- 
tigated so far. On the other hand, the investigation of 
energetical instability in the linear regime of very weak 
defects, might require times unaccessible experimentally. 

Bogoliubov Analysis. We describe the BEC dynamics 
by an effective one-dimensional GPE 
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with g = 4:Tth 2 a/m, m the atomic mass, a the s-wave 
interatomic scattering length, n the 3D density and L 
the length of the system. Moreover Vd e f (x) is the defect 
istantaneously raised at t = and O(t) the Heaviside 
theta function. Thc GPE can be linearized by writing the 
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condensate wave-function, normalised to 1 in the interval 
L, as 



if}(x,t) = -Le^-"*)/* x (2) 



being /i = p 2 /2m + gn the chemical potential. The coeffi- 
cients c q (t) are time-dependent quasi-particle amplitudes 
and U q , V g , normalized to \U q \ 2 — \V q \ 2 = 1, satisfy the 
Bogoliubov equations with dispersion 
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The momentum q represents the momentum of the quasi- 
particle excitation with respect to the momentum p of 
the condensate. The sound velocity in th e condensate 
at rest (p — 0) is given by c = -J ' gn/m. When the 
condensate momentum p < ps = mc, oj q ^ for all q, 
while for p > pb there exists a value of the momentum q 
for which oj q = 0. The time evolution of the coefficients 
c q (t) can be calculated analytically following the quasi- 
particle expansion approach developed in . We obtain 
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being I q — j- j dxe~ lqx / h Vd e f{x) the spatial Fourier 
transform of the defect. According to Eq.QJ, when 
uj q 7^ 0, the quasi-particles modes are stable: their am- 
plitudes oscillate in time remaining finite. Instead, the 
amplitude of the unstable mode with uj q — explodes, 
leading to a strong dissipation of the condensate kinetic 
energy. 

It is worth emphasizing here the main physical dif- 
ferences between energetic and dynamical instabilities. 
The latter is characterized by a complex dispersion re- 
lation so that, even in absence of perturbing potentials, 
small fluctuations in the population of the quasi-particle 
modes increase exponentially fast. On the contrary, when 
the system is energetically unstable, fluctuations remain 
small unless the condensate is perturbed by an external 
defect. In this case, Eq.QJ predicts an increase of the 
instable mode quadratically in time. 

We notice that in the Bogoliubov linear approach, the 
Landau critical velocity is indipendent of the shape and 
the strength of the defect. With a <5-like defect Vdef — 
VqS(x/L), I q = Vo for all q, so that all quasi-particle 
modes are excited. Instead, with a periodic potential 
Vdef = Vo cos(ax/h), I q = (Vo/2)S q= ± a and it is possible 
to selectively excite the two modes at q = ±a [17]]. In 
the Bogoliubov framework it is possible to predict some 
simple typical behaviours: 



1) for p ^ the spectrum becomes asymmetric and the 
dynamics, involving the two modes at q — ±a, will be 
dominated by the the quasi-particle mode with smaller 
frequency, since its amplitude is enhanced by lu q at the 
denominator; 

2) this quasi-particle excitation contains two momen- 
tum components at ±g, the relative weight of which de- 
pends on a: for small a the phonon character of the 
excitation will be strong (U q ~ V q ) and both momentum 
components become important; for large a the excitation 
is almost free-particle like and it mostly contains only one 
momentum component at q. 

From the above considerations, one can deduce that for 
large p and large a only a single momentum component 
other than the condensate will be considerably excited, 
permitting to consider a two-mode approximation, as we 
will discuss in the following. 

To inquiry about the energetical stability of the system 
using a spatially periodic defect, it is necessary to scan 
the wavevector of the periodic potential a over all possi- 
ble values. The Landau criterion can be then stated in 
the following terms: below the critical velocity, the deple- 
tion from the condensate oscillates with finite amplitude 
for all a; above the critical velocity, there exist a specific 
value of a at which the depletion grows quadratically in 
time. It is interesting to point out that a (5-like defect 
exciting simultancaously all momenta plays the role of 
the scan over all possible a mentioned just above. How- 
ever, this is of course true only in the linear regime. In 
the non linear regime mode-coupling makes things much 
more complicated and of difficult interpretation, so that 
we restrict ourselves to consider in detail only the case of 
a spatially periodic defect. 

Numerical Simulations. Bogoliubov's predictions can 
be verified numerically by integrating the full GPE. The 
numerical solutions of the GPE are well approximated by 
Bogoliubov theory until the amplitude of the depletion is 
of the order of few percents of the condensate. At longer 
times, the two evolutions differ dramatically. First of all, 
as expected, also the population of the unstable mode 
oscillates, since it cannot grow without limit. However 
stable modes might oscillate with larger amplitude than 
the unstable one, and in general they oscillate with pe- 
riods and amplitudes different from those predicted by 
Bogoliubov and strongly depending on the value of Vq. 
The main consequence is that the clear signature of the 
onset of energetic instability present in Bogoliubov the- 
ory does not exists: for each value of the momentum of 
the condensate p and each value of the wavevector a of 
the perturbing potential, the excited modes oscillate in 
time. It is not possible to claim whether for a given p 
and a given a the system is stable or unstable, but a 
comparison of the results at different a is needed. 

In order to identify a general signature of the energetic 
instability, we performed several numerical simulations 
fixing the momentum of the condensate p and changing 
the wavevector a and the strength Vq of the defect. In 
Fig. ^ we report the depletion from the condensate and 
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the population of the momentum components at q = ±a 
as a function of a for different values of p. Comparing 
the results relative to different momenta p, we observe 
different situations. For large p, the depletion presents 
a discontinuity at a certain value of a. Below a critical 
value p c this discontinuity disappears. We define oper- 
ationally the critical velocity for energetic instability as 
the smallest velocity of the condensate at which we ob- 
serve such a discontinuity in the depletion. The critical 
value strongly depends on the strength of the defect Vb, 
and, decresing Vq, the depletion asymptotically repro- 
duces the Bogoliubov prediction. 
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FIG. 1: Depletion and maximum amplitude of oscilla- 
tion of the momentum component at q = ±a for p — 
0.36,0.48,0.6,0.72 p B and Vb = 0.2 gn. The thick line is 
the depletion from the condensate, while the normal and thin 
lines below are the contributions to it from the momentum 
component at q — —a and q = +a respectively. 

At this point we can remark the relation with the previ- 
ous formulation of the Landau criterion: below the critical 
velocity, the depletion varies smoothly over all a; above 
the critical velocity, there exist a specific value of a at 
which the depletion makes a jump, which we take as the 
signature of the energetic instability. In the limit Vq — ► 0, 
the region around the discontinuity becomes sharper and 
asimptotically approaches the Bogoliubov prediction. 

From the experimental point of view, it would be eas- 
ier to study the response of the system at a fixed a as 
a function of the condensate momentum p as shown in 
Fig. |21 The jump in the depletion occurs at a momentum 
larger than the critical one, and it is a signature of the 
energetic instability of the system at such value of p. 

Following our definition, the critical velocity can be 
studied as a function of Vb, as shown in Fig. [31 For small 
Vb, we recover the Landau critical velocity, indicated by 
the dotted line at v = c. The limit of very small defect 
(Vb < 10~ 3 gn) can not be investigated numerically due 
to the very long time scale of the dynamics. The criti- 
cal velocity decreases with increasing Vq and practically 
vanishes around Vb = gn/2. For Vo > gn/2, being the 
strength of the defect larger than the interaction energy, 




FIG. 2: Depletion and maximum amplitude of oscillation of 
the momentum component at q — ±a as a function of p for 
a = 0.21,0.35,0.5 p B and V = 0.2 gn. The thick line is 
the depletion from the condensate, while the normal and thin 
lines below are the contributions to it from the momentum 
component at q — —a and q = +a respectively. 



the dynamics involves a huge number of excitations, and 
it is difficult to establish numerically the presence or the 
absence of the jump. 




FIG. 3: Critical velocity as a function of the strength of the 
defect as defined in the text. The blue dots are the values 
obtained by GP simulations, the diamonds are obtained from 
the two-mode effective potential and the dotted line indicates 
the sound velocity c. 

Nonlinear two-mode model: As we have briefly men- 
tioned above, there is a regime (large p, large a) in which 
a part from the condensate only one single momentum 
component is considerably populated and one can de- 
velop a two-mode model. We write the order parameter 
as 
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with \ak(t)\ 2 = 1. The initial conditions are dk= p = 1 
and cik^p = 0. Replacing J3J in the GPE, we ob- 
tain a set of coupled equations for the amplitudes dk(t), 
the coupling between different modes being provided by 
the presence of the defect and being proportional to its 
strength Vq. The system of equations shows that the de- 
fect couples the condensate to the modes with momen- 
tum k = pia, and to higher harmonics, with momentum 
k = p±£a, with I arbitrary integer. A drastic simplifica- 
tion in this set of coupled equations is obtained noticing 
that, when both the momentum of the condensate p and 
a are large, only the momentum k = p — a (for p > 0) 
couples significantly to the condensate. This is consis- 
tent to the picture above, which predicted that, for p 
and a large, the Bogoliubov mode with q — —a lying 
in the free-particle region of the spectrum (U q ~ 1 and 
V q ~ 0) dominates the dynamics and gives rise to a single 
momentum component at p — a. 

The equation of motion for the relative population of 
the two modes z = (N p - N p - a )/N € [—1,1] is 



becomes negligible. Based on this result, we can use the 




FIG. 4: Comparison between the results of the two-mode 
model (green lines) with the results of the GP equation (blue 
lines with dots). 



z{tf + W[z(t)} = const., (6) 

with 

w[z) = -^o 2 (l - * 2 ) + [A£(l - z) gn(l z 2 )/2] 2 (?) 

where AE — p 2 /2m — (p — a) 2 /2m, and the initial con- 
dition are zo = z(0) = 1 and i(0) = 0. These equations 
are similar to the Josephson equations for the double 
well and many other two- mode nonlinear systems |19| . 
With a mechanical analogy, the dynamics of the system 
is described by the motion of a particle in the potential 
W(z) in Eq. |JBJ. When the condensate has momentum 
p smaller than the critical value, the effective potential 
displays a double- well shape for every value of a. The 
potential barrier remains higher than the initial energy 
of the system so that the motion is always localized in 
the right well. On the other hand if the momentum p 
of the condensate is larger than the critical value, by de- 
creasing a the potential W(z) displays a transition from 
a double-well to a single-well, i.e. the potential barrier 
becomes smaller than the initial energy of the system. In 
this case, the population of the excited mode suddenly 
increases, qualitatively reproducing the jumps obtained 
with the GP simulation shown in Figs. l|ll2|l . 

We have checked the capability of the two-mode model 
of reproducing the results of the GP simulations, and 
confirmed that it works well for large p and large a, as 
shown in Fig. 

Close to the critical velocity, the two-mode model does 
not work quantitatively well and overestimates the crit- 
ical velocity, as shown in Fig. [3] However it reproduces 
the decrease of the critical velocity as a function of Vq 
and gives correctly the point at which the critical velocity 



two-mode model to calculate the critical defect strength, 
at which the critical velocity vanishes. To that aim, we 
consider small p and small a, i.e. small AE. Within this 
assumption the effective potential can be expanded and 
becomes 



W(z) « ~V 2 (l-z 2 ) + (gn/2) 2 (l-z 2 ) 2 + (8) 
-gnAE(l - z 2 )(l - z). 

When Vo — gn/2, the effective potential undergoes a 
double to single-well transition for each value of p: for 
a > 2p, AE < and we have that W(0) > W(l); for 
a e [0, 2p], AE > and we obtain that W(0) < W(l). 
This means that at a ~ 2p, W(z) undergoes a transition 
from double to single-well, and hence the population of 
the excited mode undergoes a jump, signature of the en- 
ergetic instability. 

Conclusions. The Landau criterion for superfluidity 
predicts that a weakly interacting Bose-Einstein conden- 
sate can flow without dissipation when its group velocity 
is smaller than the sound velocity. We found that the 
Landau criterion remains valid only in the linear (Bo- 
goliubov) regime or, in other words, for vanishingly small 
perturbations. Indeed, the superfluidity critical velocity 
rapidly decreases while increasing the strength of the de- 
fect, up to a critical value of the perturbation where the 
flow is always dissipative. This finding can be tested ex- 
perimentally with spatially periodic perturbations, and 
can qualitatively explain the critical velocities lower than 
sound observed with condensates flowing against large 
defects. 
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